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Abstract 

We examine the effects of new physics on the Higgs sector of the standard model, focusing on 
the effects on the Higgs self couplings. We demonstrate that a low mass higgs, < 2 nit, can 
have a strong effective self coupling due to the effects of a new interaction at a TeV. We investigate 
the possibility that the first evidence of such an interaction could be a higgs-higgs bound state. 
To this end, we construct an effective field theory formalism to examine the physics of such a low 
mass higgs boson. We explore the possibility of a non relativistic bound state of the higgs field 
(Higgsium) at LHC and construct a non relativistic effective field theory of the higgs sector that 
is appropriate for such studies (NRHET). 



I. INTRODUCTION 



Currently, global fits to all precision electro weak give the higgs mass to be 113+^GeV 
with an upper bound given by < 241 GeV at 95% CL (see, e.g., J. Erler and P. Langacker 
in Sec. 10 of Ref. [l]). LEP has also placed a lower bound limit of > 114.4 GeV ^J. 
Assuming the Standard Model (SM) of electroweak interactions, one expects that the higgs 
will soon be found at LHC. 

However, there are at least two reasons why the SM with a single higgs doublet is ex- 
pected to be incomplete. The first is the triviality problem. This asserts that the higgs self 
interaction, and hence its mass, must vanish unless the theory has a finite cut-off. Triviality 
has been rigorously established only for simpler models, but it is widely believed to hold for 
the SM higgs. The other is the hierarchy problem that quadratic divergences need to be 
finely tuned to keep the scale of electroweak breaking smaller than the natural cut-off of the 
theory (which in the absence of new physics would be the Planck scale). 

In this paper, we investigate the effects that new physics, invoked to cure these problems, 
may have on the higgs sector of the SM. We assume that the scale of the masses of new 
quanta, Ai, is sufficiently higher than the scale of electro- weak symmetry breaking (v ~ 
246 GeV) so that the quanta of the unknown new physics can be integrated out. As we want 
this new physics to address the hierarchy and triviality problems, and for phenomenological 
reasons, we are interested in new physics where M. ~ TeV. The resulting low energy 
effective theory is the one higgs standard model supplemented with non-renormalizable local 
operators, of dimension D > 4, which are constructed of standard model fields invariant 
under the SU(3) x SU(2) x U(l) gauge symmetry. This approach has been applied to 
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and has recently been the subject of further 
15j . The advantage of this approach is that it is model 



independent: any new physics scheme that results in a low energy spectrum coinciding 
with the SM's can be described in this way. The disadvantage is that the new physics is 
parametrized in terms of several arbitrary parameters, the coefficients of higher dimension 
operators, and nothing is known a priori about these coefficients. 

For a particular extension to the standard model, consistency requires that fits such as [l| 
be reconsidered with the new operators, severely relaxing the constraint on the higgs m ass. 



In fact, it has been shown that the effect of higher dimension operators 16, 
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elminate the mass limit on the higgs. While more exotic possibilities are tantalizing, in this 
paper, we focus on the possibility that the new physics integrated out is strongly interacting 
and effecting the higgs sector above the scale M. while the higgs itself has a relatively low 
mass rrifi < 2m t . 

Various bounds can be placed on M. from low energy experiments. In particu- 
lar, flavor changing neutral current bounds such as those arising from K° — K° mix- 
ing impose strong constraints, M. > 10 4 TeV. These bounds can be relaxed by re- 
stricting the higher dimensional operator basis through adopting the MFV hypothesis 
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281 ] . This allows one to consider M. to be a few TeV 



while naturally suppressing FCNC. 

However, even utilizing the MFV hypothesis to justify new physics at a TeV, higher 
dimensional corrections to the standard model could exist that modify the relation mw — 
mz cos^iy. This relationship is experimentally required to be respected to a fraction of a 
percent. The PDG quotes po = 1-0002+q;qqq4 for the global fit 29J of precision electro-weak 
observables. This fact motivates the consideration of new physics being integrated out that 
preserves po ~ 1 naturally, even with possible strong dynamics effecting the higgs at the 
scale M.. 

This can be accomplished assuming an approximate custodial SU(2)c symmetry 
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321 ] . where the weak SU(2) gauge vector bosons transform as a triplet and the higgs field 



transforms as a triplet and a singlet. The higgs vacum expectation value is in the singlet 
representation of SU(2)c , so the approximate symmetry is explicitly realized, and is ex- 
plicitly broken only by isospin splitting of fermion yukawa couplings and by hypercharge. 
We require that the operator extensions to Csm respect this SU(2) C symmetry up to hyper- 
charge and Yukawa coupling violations, as in the standard model. Operators that break the 
custodial symmetry are allowed but their coefficients are taken to be naturally suppressed. 

In the SM, the higgs cubic and quartic couplings are not independent parameters, but 
given in terms of the higgs vacuum expectation value v and mass rrih- The obvious immediate 
effect of D > 4 operators is to shift all these quantities in independent ways, so that 
effectively the higgs cubic and quartic couplings become independent parameters. Of course, 
the shift from the SM values is somewhat restricted, of order (v 2 /Ai 2 )C, where C is a 
dimensionless coefficient, C ~ 1. The effect on single higgs production rates of modifications 
to the coupling of the higgs to weak vector bosons or to itself was investigated in Ref. 11]. 
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The modification of higgs decay widths and this general class of models was also examined 



in 



33]. 

In this paper we address the question of whether a bound state of higgs particles can 



form. In the SM a higgs bound state forms only if the higgs is very heavy [341. |35|. There is a 
competition between the repulsive interaction of the quartic coupling, Ai, and the attractive 
interaction of the higgs exchange (between higgs particles) which is determined by the cubic 
coupling, Aii>. For large enough coupling the exchange interaction is strong enough to 
produce binding, but since the mass, yX\v, is a l so given in terms of the coupling, the higgs 
mass is large. In the effective theory context the three parameters (mass and cubic and 
quartic couplings), are independent and a bound state is possible for smaller higgs mass. 
The question of detail becomes, how is the bound on the higgs mass for a bound state to 
form relaxed by the coefficients of D = 6 operators? Can one have a bound state of light 
higgses? We find that the effect of these operators can be significant, allowing for bound 
states for much lighter higgs particles. Discovery of such bound states would give valuable 
information on the scale of new physics. 

There is no know solution to the bound state problem for identical scalar particles inter- 
acting via cubic and quartic interactions. The higgs bound state problem has been addressed 
using different approximations, the N/ D method is Ref. 



homogeneous Bethe-Salpeter equation in Ref. 



341 ] and a truncated version of the 



351 ] . Our aim here is to find a necessary con- 



dition on the coupling for which a non relativistic (NR) bound state may form. To this end 
we introduce a new method. We propose to study the formation of the bound state in a 
non-relativistic effective theory for higgs-higgs interactions. 

We begin by listing the D = 6 operators of the effective theory. We take two approaches. 
In the first, linear realization, we consider operators that can be built out of the higgs 
doublet and the fields in the gauge sector of the SM. Our primary interest here is in the 
higgs sector per se, so we focus on higgs self interactions. The second approach, uses a 
non-linear realization of the symmetry. Since the higgs field is intimately connected to the 
symmetry breaking of the SM gauge symmetry, it is natural to expect that below the scale of 
new physics the effects of symmetry breaking are already apparent. Were the higgs mass as 
large as the scale of new physics, the SM would be supplemented not with a higgs doublet 
but with a triplet of would-be goldstone bosons that are eaten by the W and Z vector 
bosons. The higgs, if somewhat lighter than the scale Ai, would appear as a singlet under 
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the gauge symmetry. 

We then proceed to construct the effective theory at low energies. If mh ~ 2m t) one can 
incorporate the virtual effects of the top by integrating it out and constructing a top-less 
effective theory. In order to investigate the minimal coupling for which a NR higgs-higgs 
bound state may form we then construct a non-relativistic higgs effective theory, and proceed 
to determine this condition. 



II. HIGGS EFFECTIVE FIELD THEORY: LINEAR REALIZATION 
A. The D = 6 Custodial SU(2) Higgs Sector 

The Lagrangian density of the standard model containing the higgs field 1 is given by 

Cl = (D^y (D.^-V^) (1) 

where <fi is the higgs scalar doublet. The covariant derivative of the <fi field is given by 

D» = ld,-i 9 -±B,-ig 2 ^Wl (2) 

where a 1 are the pauli matrices, W:, are the SU(2) and U(l) SM gauge bosons and the 
hypercharge of 1/2 has been assigned to the higgs. The higgs potential at tree level is given 
by 

V(0) = -m 2 t 0+y (0V) 2 - (3) 

No dimension five operator can be constructed out of higgs fields and covariant derivatives 
that satisfies Lorentz symmetry and the standard model's gauge symmetry. 2 Utilizing the 
equation of motion of the higgs field and partial integration the number of dimension six 
operators is reduced. The effective Lagrangian density of the extended standard model is 
given by 



1 We have omitted Yukawa interactions with fcrmions here. 

2 To satisfy Lorentz invariance an even number of covariant derivatives are required. To be invariant under 
the SU(2) x U(l) gauge group the operator must be bilinear in 0t and <fi. 
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where the dimension six operators that preserve the symmetries of the standard model and 
custodial SU(2)c in the Higgs sector are given by 

4 = C<^ (0f (0t 0) + C<J (0t0) (Z^t (D^)-^ (0t0) 3 . ( 5 ) 
Note that the operators considered here preserve custodial symmetry and can result 



from tree level topologies in the underlying theory. [36|] As such, these operators need not be 
suppressed by loop factors of 1/167T 2 or proportional to a small custodial symmetry breaking 
parameter. For these reasons these operators are expected to have the dominant effects on 
the higgs self couplings and we take their coefficients to be 0(1). There is only one operator 
in the Higgs sector that violates custodial symmetry and could come from an underlying tree 
topology, {(jr D^cf)) 2 . The underlying topology in this case determines that the symmetry 



breaking parameter is given by g\. The coefficient of this operator has been determined 37] 
to be C < 4 x 10 -3 where we have used A = 1 TeV. We neglect this operator. 

We expand the higgs field about its vacuum expectation value with (h(x)) = and treat 
v 2 1 M 2 as a small perturbation. We expand the field as usual around a vacuum expectation 
value v so that 

^ (.;*,)■ 

Here U(x) = e l ^ a ^ C7a ^ v and the would-be goldstone boson fields of the broken symmetry are 
£ a . In unitary gauge, the gauge transformation can be used to remove the goldstone boson 
fields. We then redefine the higgs field (h) so that the kinetic term is normalized to 1/2, 
using the field redefinition 

h h ' m 

(l + 2Cf)V2' V) 

where Off = (v 2 /A4 2 )(C^ + \C 2 ). The effective Lagrangian density is given, in terms of the 
rescaled field, by 

C% + -j± = -d> l tid li ti- V eff (ti) + c$ o % $ + c ww o ww + c zz o zz 

~T ^h'WW w h'WW ' Z Z ^h' Z Zi \° ) 
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summed over i,j such that i + j = 2 , where 
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The effective potential is 
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which is written in terms of the rescaled mass term and the effective couplings, which are 
given by 
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2 """A* 2 ' ~ V A4 4/ ' 

We will suppress the prime superscript on the higgs field for the remainder of the paper. 



(12) 
(13) 
(14) 



B. D = 6 SM Field Strength Operators 



The operators that can be constructed out of the higgs scalar doublet and the field 
strengths (or duals) of the standard model are as follows. We restrict our attention to those 
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operators listed in p, ll5[ that preserve the SU(2)c custodial symmetry: 

Here G^, W 7 "^ and 5^ stand for the field strength tensors of the SU(3) x SU(2) x Z7(l) 
gauge bosons, and a tilde denotes the dual field strengths, F fll/ = e^ u Acr F Xa /2. Note that 
the operator that is proportional to the 5" parameter given by 

-^(^V^)r% (16) 

violates custodial symmetry and is naturally suppressed in our approach 3 . 
C. D = 6 Fermion Sector 

Operators of dimension 5 and higher that couple the higgs to fermions, or purely fermionic 
operators, can give rise to unacceptably large flavor changing neutral currents (FCNC). If 
the coefficient of such operators are generically of order 1 the scale of new physics must be 
taken to be Ai ^ 10 4 TeV in order to suppress FCNC effects. We adopt the Minimal Flavor 
Violation hypothesis (MFV) to naturally suppress the dangerous operators while maintaining 
a low scale of new physics, M. > 1 TeV. In the absence of quark and lepton masses the SM 
has a large flavor symmetry group, Gp = SU(3) 5 . The MFV asserts that there is a unique 
source of breaking of this symmetry. All operators that break the symmetry must transform 
precisely the same way under Gf- As a result FCNC operators are suppressed by the 
familiar factors of the Kobayashi-Maskawa (CKM) matrix in the quark sector and by the 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix and small neutrino masses in the lepton 
sector. 

Since the effects of fermionic operators are not needed for the rest of this investigation, 
we do not list the operators. The interested reader can find a complete description of the 
operators and their effects in 6|]. 



See Appendix A 
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III. HIGGS EFFECTIVE FIELD THEORY: NON LINEAR REALIZATION 



The construction in the previous section assumes that the field content of the effective 
theory includes a higgs doublet. This is not necessary. If the electroweak symmetry is 
spontaneously broken by a strong interaction the spectrum below the scale of this new 
physics does not have to be described by a higgs doublet field, beyond the SM fields. Only 
fields describing the would-be goldstone bosons need be introduced. Such higgs-less theories 

]. However, if the higgs particle is somewhat lighter 



have been discussed in the literature 



than the scale of new physics it has to be incorporated in the low energy description and 
symmetry alone does not dictate that it appears as a member of an iso-doublet. It is sufficient 
to have the goldstone bosons realize the broken symmetry non-linearly, and the higgs field 
is then a singlet under the symmetry. 

The situation is entirely analogous to the case of 7r's and the a in QCD. A phenomeno- 
logical Lagrangian density describing it and a interactions does not have to be a linear 
realization of the chiral SU(2) x SU(2) symmetry. Instead, the 7r-fields have a better de- 
scription through a non-linear chiral Lagrangian. Then the a can be included through 
interactions that satisfy the non-linearly realized symmetry and the usual rules for naive 



dimensional analysis 39]. 



In the non-linear realization, the Lagrangian density in Eq. (TjQ) is replaced by 

£ NL = yTrD^D^V + \d^h - V(h), (17) 

where the would-be goldstone bosons £ a appear through the matrix U(x) = e l ^ a<yX ^ <Ta / v that 
transforms under SU(2)l x SU(2)r linearly, U — > LXJBJ , and h is a singlet field, describing 
the higgs particle. Custodial symmetry SU(2)c is the diagonal subgroup of SU(2)l x SU(2)r 
and the higgs field is invariant under it. 4 

This Lagrangian is supplemented by higher order terms suppressed by powers of Ai. In 
the case of the higgs potential, this can be included simply as 

V(h) = M 4 f(h/M), (18) 

where f(x) is an arbitrary function with a minimum at zero. The mass and couplings of the 



The custodial symmetry is discussed in more detail in Appendix A. 
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higgs are given in terms of this dimensionless function by 

ml = M 2 f("\0), (19) 
vX* f = MfW(0), (20) 
Af = /W(0). (21) 

It is not a surprise that in the non-linear realization of the symmetry the couplings and mass 
are completely independent, and that they are all naturally of order 1 times the appropriate 
power of the dimensionfull scale, Ai. The natural scale for the higgs mass is At, and we 
are considering here the class of theories for which f^"\0) happens to be small, while higher 
derivatives may remain of order 1. We stress that the natural scale for the cubic coupling is 
Ai. Unless the mechanism (or numerical accident) that keeps the higgs mass small compared 
to Ai also acts to suppress the cubic coupling, one must naturally expect X e J^ ~ Ai/v ^> 1. 
We will also need the corrections to the derivative interactions. We write, generally, 



eff i e ff 

Ci- + c 2 JJ — 
V v z 



d»hd,h-± m lh*- V -^h*-^h* + ... (22) 



In the linear realization the derivative interaction couplings are related, = \c 2 = 

Off = {v 2 1 ' Ai 2 ){C^ + C|/4), but in the non-linear realization they are independent. And, 
as in the case with A3 naive dimensional scaling gives an enhancement of that could 
arise from the non-perturbative dynamics of the symmetry breaking sector. Naively, ~ 
(v/At), which is enhanced over the linear realization value by a power of (Ai/v). 

As we mentioned earlier, non-linear realizations have been extensively studied for higgs- 
less theories, but have been neglected in studies including a light higgs. There are two 
important consequences of the non-linear realization outside the pure higgs sector that we 
point out here. It has been noted that significant corrections to the coupling of a higgs to 
gluons are possible from D > 4 operators. The modifications can be large because there 
is no SM contribution at tree level. In the linear realization there is a D = 6 operator 
that contributes at tree level, and therefore competes with the SM one loop, top mediated 
amplitude: 

<G 'uvG [ uv, (23) 



Note that the linear realization implies a relation between the one and two higgs couplings 
to two gluons. However, in the non-linear realization the two couplings are completely 
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independent 

h 



Cl M +C2 M-i) Ga ^ (24) 
In Ref. |40j it was noted that a heavy quark with Yukawa coupling A — > oo produces a 

coupling of two gluons to one or more higgs particles that cannot be described by the 

effective theory operator in ( |23|) . Instead a nonpolinomial interaction was introduced to 

describe this effect, 



There is no problem accommodating such interactions in the non-linear realization, by 



g ln(l + h/v)G%G% = %L[h/v + (h/vf + ■■■ ]G%G« U . (25) 



In much of what follows we implicitly assume the linear realization. However, results in 
terms of the arbitrary parameters m^, A3 and A4 can be interpreted readily as arising 
from the non-linear realization. 

IV. A LOW ENERGY EFFECTIVE THEORY FOR THE HIGGS 

In this section we will construct an effective theory for the light higgs, integrating out 
momentum modes heavier than the higgs. This is useful in discussing physical effects with 
a typical energy of order of the higgs mass. In particular, we integrate out the top quark. 
As the coupling of the top quark to the higgs is fairly large, we would like to estimate the 
effects of the top quark on the possibility of forming a higgs bound state. If the top quark 
mass is much heavier than the higgs it is appropriate and convenient to describe the higgs 
self interactions in a top-less theory. When the top quark has been integrated out, its effects 
are accounted for through modifications of coupling constants and mass of the higgs. 

While this is clearly appropriate when the top quark is much larger than the higgs mass, 
we use this approximation even when the higgs is slightly heavier than the top. For < 
2m t the approximation is known, ipso facto, to work better than one would expect. This is 
due, in part, to the fact that there is no non- analytic dependence on the mass since the higgs 



is the pseudo-goldstone boson of spontaneously broken scale invariance|4ll. l42l. l43j] . It is also 
known that soft gluon effects are large and correctly reproduced by the effective theory 44]. 

For single and double higgs production, comparisons between the the full theory calcu- 
lation and the effective top-less theory find that the latter is a good approximation for the 
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total rate for rrih < 2m t . For example, with the appropriate K factor, the resulting topless 



he full NLO 



44]. 



effective field theory calculated to two loops is known to accurately describe t 
result for g g — > h to better than 5% accuracy in the full range < < 2 m t 

As another concrete example, consider the higgs mass dependence in the higgs IPI self- 
energy. The first graph of Fig. [T] is the contribution of the top quark to the IPI self-energy 
which we label — iU(p 2 ). At 1-loop we find 



N c m\ 
An 2 v 2 



dx ( 1 



3 log 



li- 



ra 2 — x(l — x)p' 2 



(26) 



where we have used the MS subtraction scheme. The quantity 1 — U(0) /U(p 2 ) (at fi = m t ) 
never exceeds 30% when ^/pP■ ranges from zero to 2m t . 

For these reasons we consider it appropriate to integrate out the top quark for < 2m t 
in this initial study. When >> m t , these corrections should be taken only as indications 
of the size of virtual top effects. While the approximation of neglecting higher order terms 
in the p 2 /m 2 expansion is known to work better than expected for the applications we will 
consider, there is no guarantee that it will work well for processes not considered here. 
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A. Running to mt 

The coefficients of the D = 6 operators at the scale A4 are unknown. We are assuming 
that the new physics couples to the higgs field and is strongly interacting at the scale M.. 
In this context, it is natural to take 

C;(M),\ 2 (M)~l. (27) 

Similarly it is natural to assume that the coefficients of the D > 4 operators that couple the 
higgs to other fields, like those in Eq. (fl~5]) or those that couple the higgs to quarks while 
satisfying the MFV hypothesis, are all order unity. 

The anomalous dimensions of the extended operator basis can be determined system- 
atically. This is beyond the scope of this paper. But the effect of the running is easy to 
understand. With minimal subtraction the calculation of the running of coefficients of higher 
dimension operators can be done in the symmetric, massless phase. There is operator mixing 
among the D = 6 operators with common quantum numbers. The anomalous dimension 
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matrix is a function of the relevant couplings (Ai, g±, gi and the top quark Yukawa, At). 
The running is always proportional to these coefficients so the effect is roughly of the form 

C\{m t )~C\{M) (l + ^bg^)), (28) 

where mixing is implicit, and c\a stands for a linear combination of Ai and the squares of 
gx, g 2 and A t . 

Since log (m t /Ai) ~ 1 and the coefficients c\ ~ 1 the running produces a small, calculable 
shift in the unknown coefficients. Hence, we continue to take the unknown Wilson coefficients 
at the scale m t to be ~ 1 . 

At m t the top quark is integrated out and this produces a different effect, a shift in the 
C^(m t ) by a C^- independent amount. This can be numerically significant, and we estimate 
this next. Note that once the top is integrated out we continue to run down to the mass of 
the higgs scalar nih- The effect of the running of these coefficients from m t to mh is again 
small, so we take 

C\ (m h ) , A 2 (m h ) ~ 1. (29) 



B. Integrating out the top quark 

Integrating out the top leads to further corrections to the higgs sector of the standard 
model. The top mass is a result of symmetry breaking, so the resulting effective theory is 
better presented in unitary gauge, as in (jSJ) and (iTTj) . In unitary gauge, the top mass term 
and coupling to the higgs is given by 

C Y = -m t q t q t M + ^ . (30) 

We begin by considering effects on the higgs self- couplings. Figure [1] show the Feynman 
graphs that contribute to modifications of the higgs self-couplings. The solid line denotes a 
top quark, the dashed external lines denote the higgs. 

We perform the calculation to lowest order in p 2 /mf. Some details of the computation 
are given in the appendix. The effect of these corrections is to further modify the effective 
potential of the higgs scalar field h. The effective couplings and mass term of Eqs. (fl2"]) - (fH|) 
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FIG. 1: Integrating out the top quark, 
are shifted by these corrections, and are now are given by 




As emphasized above, these corrections are not multiplicative, that is, they are present 
even for A2 = = 0. Whether they are important depends on the scale and strength of 
the new physics. The condition 

< 34 » 

is satisfied for A2 ~ 1 when M. ~ 2irv = 1.6 TeV. So the corrections are numerically 
comparable to these new physics terms. Similarly, for Ai ~ 1 the condition 

Cf*t = ^(C5 + JCj)*'~(^)pf < 35 > 
still requires M w 1.6 TeV for C\ + \C\ ~ 1. 
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C. Corrections to Field Strength operators 



Integrating out the top quark also results in effective operators of the higgs field and 
the SM field strengths. The dominant SM production mechanisms for the higgs at LHC is 
the gluon fusion process g g — > h. We restrict our attention to such operators that effect 
the production processes of the higgs through gluon fusion. Figure [2] shows the 1-loop 
Feynman diagram for the top contribution to gg — > h. For a higgs with m/, < 2m t , the 
expected production cross section of the g g — > h process has been determined up to NNLO 
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52j. For SM gluon fusion, the single higgs production mechanism is given by the 



m t — > oo effective Lagrangian density comprised of a dimension five operator 

, , h 
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CGGh ( a s 



V 



(36) 



where t 
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re coefficient is given in the MS scheme, in terms of a s for five active flavors, by 
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(37) 




FIG. 2: The gluon fusion g g — > h production process. The production process through the effective 
local operators in shown in the second column. The effective local operators come from integrating 
out the top quark and new physics at M. 

The interactions in the effective Lagrangian of Eq. (fTBT) also contribute to single higgs pro- 
duction through gluon fusion. Combining results, at the scale m^, the effective Lagrangian 
density for single higgs production is given by 

= Cf Gh \ G% G AfiU + Cf Gh \ G% G AtJ,v , (38) 

where 

C e JL = C GGh -2na s c G ^ (39) 

Cf Gh = -2na s c G -^. (40) 
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Assuming that the new physics degrees of freedom carry the SU(3) gauge charge, the 
Wilson Coefficients cq, cq will be approximately the same size as the coefficients C^, A2 we 
are interested in. If the new physics degrees of freedom are charged under SU(2) x U(l) but 
not SU(3), below the scales A4, m t effective local operators of this form will still be induced. 
However, the corresponding Wilson Coefficients will be suppressed by factors of 16 tt 2 . 

The effect of these interactions on higgs production rates was examined in [l4]. Note that 
in the standard model, contributions to the operator G A v G Ailu are highly suppressed j|39] 
and therefore neglected. 

^umrouuuuwx — ' — 7 




1 c, o, 



FIG. 3: The gluon fusion g g — > hh production process and the effective local operators. 



The production process of two higgs in the standard model is shown in Fig. [3J In analogy 
with the single higgs production case we characterize the process in the effective theory by 



h 2 h 2 

^ GGhh „,2 /if" ' GGhh „,2 M v " 5 



where the coefficients are given by 



r< e ff — n 1 

^GGhh — ^GGhh 



c G 7ia s v' 



M- 



^GGhh 



C G TT a s V z 
M 2 



Here the top quark contribution is 56| 

^GGhhif^ 2 ) = 



(41) 

(42) 
(43) 

(44) 



12 7T 48 TT 2 

The Wilson coefficients for two higgs production in the effective theory is not suppressed 
relative to the corresponding Wilson coefficient for single higgs production. Note that unlike 
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the case of single higgs production the expansion in p 2 jm\ does not, in general, have kine- 
matics such that p 2 /m 2 ~ m 2 h /m 2 . In two higgs production, higher order terms in p 2 jm\ 
have p 2 = s,t,u and in general (s,t,u)/m 2 is not small. We calculate the next order in 
the expansion of p 2 /m 2 in Appendix [B] These terms are neglected, and our application of 
the expansion is valid for finite values of m t due to our interest in establishing a necessary 
condition for a NR bound state to form. The kinematics for the production of a NR bound 
state at threshold dictate (s,t,u)/m 2 ~ m\jm\. 



V. PHENOMENOLOGY OF HIGGS EFFECTIVE THEORY 



A. The Magnitude of Self Couplings 

The effect of the D = 6 operators in the effective potential cause corrections to the 
three and four point contact interactions and m^. To illustrate that the induced effects on 
the higgs sector are under control, consider extending the effective potential with a single 
dimension eight term. We find the following while neglecting the effects of integrating out 
the top quark 

^ = Aj (l-2Cf) 2 + 3xl(T 2 A 2 (l-2Cf)+4.5xHT 4 A 3 , (45) 
\ e /f = 3Ai (l-3Cf + 7.5(Cf) 2 ) + 1.5 x KT 1 A 2 (l-3Cf) +3.1 x KT 3 A 3 , (46) 
X e / f = 3Ai (l-4Cf + 12 (Cf) 2 )+ 4.5x10-^2 (l - 4Cf ) + 1.6 x KT 2 A 3 . (47) 

From which one sees we are examining the potential of the theory in a controlled expansion, 
even for Ai ~ 1 TeV. 

Eliminating the self-coupling Ai in favor of the higgs mass, we can write for the effective 
cubic and quartic higgs-self couplings, 

m\ v 2 7N C mf 
~v 2+ 2 M 2 ~4^ 2 ~ 4 

and 



A 3 - 3 (1 - C h ) — + A 2 -— -j-5-t, (48) 



A^ = 3(l-2Cf)^ + 6A 2 ^-^4 (49) 

With v = 246 GeV, m t = 174 GeV and M = 1 TeV, and taking m h = v/2, these are 

A^- f/ = 0.62 - 0.05(C; + \C\) + 0.06A 2 , (50) 
\ff = 0.39 - 0.09(C; + \C 2 ) + 0.36A 2 . (51) 
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For negative A2 of order one one can greatly reduce the repulsive contact interaction, A4 , 
in a putative higgs-higgs bound state. Of course, this comes at the price of reducing the 
attractive interaction, governed by A3 . 



B. g g — > hh Production 



From our results in Section IIVCI the production of two higgs in our effective theory 
framework is straightforward to write down. The contributions to the amplitude are shown 
in Fig. HI 





FIG. 4: The two higgs production process in the effective theory. 



The amplitude for two higgs production, to 0(a s ), is given by 

(hh\iA\A a (P 1 )A l3 (P 2 )) = (hh\iA 1 \A a (P 1 )A /3 (P 2 )) + (hh\iA 2 \A a (P 1 )A l3 (P 2 )) (52) 



where we have 

{ t A 2 y? = 4iC P (CJ^ t jr"(Pi,-P2), 
where 

pa/3 



\e// 

V X 3 JJ 



\c* (Pi+pi+iP.+p^y 

(53) 



f ap (P u P 2 ) = Pf P l i + Pf P 2 ° - 2 g al3 P • P 2 



(54) 



Using two hi 



examined in 



ggs production as a test of the cubic self coupling of the higgs has been 



60] where testing for the MSSM with this signal was investigated. Ref 



advocated the examination of g g — > h h to compare the one and two higgs production 
coefficients in Eq. ( |24|) since the naive relation between the two coefficients could be upset 
by the presence of novel operators like that in (1251) . As we have discussed, in the non linear 
realization of broken electro- weak symmetry, the relationship between g g — >• hh production 
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and g g — > h production is not fixed as in the linear realization. Any deviation from the SM 
value for g g — > h h must be interpreted with care. The g g — > hh production rate in our 
effective theory construction (in the linear realization) depends on at least six unknowns, 



namely, Ji4, X<2,C^,C^,co,Cg- The effects of the operator advocated in [40J increase the 
number of unknown parameters still further. 

Clearly two higgs production is an important signal to test the higgs mechanism in the 
standard model. The cross section of g g — > h h is suppressed compared to the cross section 
of gg — > h by a factor of 1000, due to the effects of parton distribution functions and 
phase space suppression jsoj] . The cross section falls from 50 fb to 10 fb as the higgs ranges 
in mass from 100 GeV to 200 GeV. Thus once LHC enters its high luminosity running of 
100 fb~ 1 /Year one can expect roughly 1000 events per year. A significant excess or deficit of 
this signal should be observable. However, the reconstruction of exactly what form of new 
physics is present requires more information. 

One could obtain more information on the unknown parameters involved by further 
probes of the physics of the self interaction of the higgs. In the remainder of the paper 
we examine the sensitivity of a higgs bound state (Higgsium) in the appropriate low energy 
effective field theory on TeV scale physics to these parameters. If a bound state forms, one 
can use the properties of the bound state such as its binding energy, as an a probe of the 
physics above the scale M.. 

C. Higgsium: Production and Decay time 

To get some rough understanding of the conditions under which a higgs-higgs bound state 
may form, consider the non-relativistic Schrodinger equation 

[-V? + V(r) - E] ^(r) = 0, (55) 

with the potential from a yukawa exchange and a contact interaction, 

„2 p -m h r 

V(r) = -■%-- + K5 3 (r). (56) 

We are interested in the case g ~ A3 and k ~ A4 as a non relativistic approximation of 
the higgs self interactions. Neglect for now the contact interaction. It is well known that 
the Yukawa potential produces bound states provided 

£ > 1.7, (57) 
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Neglecting new physics effects, 



A5"*3=*, (58) 



non-relativistic bound states could be expected for 

m h > l.2v. (59) 

The effect of TeV scale new physics changes the relationship between the mass and the 
coupling. The above Yukawa bound state condition is modified to 

m h > [1.54 + 0.09(C; + |C|) - 0.02 A 2 ]^ v. (60) 

This demonstrates the point that if the higgs self coupling is significantly stronger due to 
strong TeV scale new physics that contributes large Wilson coefficients, then a low energy 
signal of this higher scale physics might be a NR bound state formed by two higgs. 

However, one can see that it is difficult to realize the NR bound state condition when 
we identify the couplings in this Schrodinger equation with our effective couplings. This 
identification is in fact incorrect. We will demonstrate in Section [VI] that the correct NR 
limit of the higgs sector is described by a Lagrangian containing only contact interactions 
and higher derivative operators. 

The formation time of the bound state can be approximated by the ratio of 4 Rq/u where 
i?o is the characteristic radius of the NR bound state and u is the relative velocity of the 



6l| 



two higgs. This is roughly the period of oscillation for S wave states 

For a NR bound state we can approximate the relative momenta of the two higgs by 
p ~ rrih u so that 

4#o 



T f ' 

u 



rrihU 2 



(61) 



The SM higgs decays predominantly via h — ► bb pairs through Yukawa interactions if 
114.4 < rrih <C 2 Mz- We take these decays as dictating the decay width of Higgsium. 
Neglecting the effects of our new operators, this decay has the decay width 

m? 3m h ( m 2 f \ 3/2 



v 



4:71 \ m h/ 
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This gives an approximate decay time 



4 7T V 2 



n = - 2- ( 63 ) 

3m h ml 

The condition that the bound state has time to form is that tj < r& which can be satisfied 
for 

«* >!&)■ <-) 

Thus a non relativistic bound state has time to form before it decays. Above 135 GeV and 
below the threshold of W + W~ production, the dominant decay of the higgs is through a 
virtual W pair, h — > WW*. Above rrih > 2mw the decay into W + W~ predominates and 
the decay width is given by 



w 



TTh 1 

-T^T" VI ~a w (4-4a w + 3a 2 w ) , (65) 
V 62, ix 



where aw = ^ m w/ m h using the notation of 29(. Comparing the formation and decay time 
for a higgs whose mass is above the threshold of W + W~ production we find 

u 2 > , 1 (66) 
rj(m h ,m w ) 8% \ v 2 J 

where r](m h ,m w ) ~ 1. 

The lower bounds on u in either case are compatible with the NR approximation for the 
full range of higgs masses we consider. A relativistic bound state is also possible in either 
case, although an approximation scheme that can estimate its formation time is lacking. In 
the remainder of the paper we focus on the possibility of a NR bound state being formed by 
a relatively light higgs, m h < 2m t , due to our treatment of the top quark. We also note that 
a NR bound state may also have observable effects on the spectrum of two higgs production 
even if a bound state does not fully form as in the case of top quark pair production near 
threshold in e + e~ collisions 61]. 



VI. NON RELATIVISTIC HIGGS EFFECTIVE THEORY 

If the two higgs are created with small relative velocity and form a non relativistic bound 
state it is appropriate to describe the physics of this state with a non relativistic effective 
field theory of the higgs sector. We refer to our effective theory derived in Section [J through 

21 



Section HV CI as higgs Effective Theory (HET) and now match onto a non-relativistic version 
of this theory (NRHET) where we take the c — > oo limit of the scalar field Lagrangian 
density of HET. Recall the Lagrangian density is of the form 



2 



tth tt h 2 

eff i eff 
C{ - + C 2 J ^7 

1> IT 



Phdrh-\mlh*-^h>-^h* + 0^). 



We wish to construct the non relativistic limit of this Lagrangian density systematically, 
retaining H — 1 and making factors of c explicit with [c] ~ The dimensioned 

quantities can be expressed in units of length [x] and time [t]. As h — 1, we still have 
[i?] ~ l/[i] and [p] ~ As the action S = J dtd 3 x£ is dimensionless, we have [C] ~ 

[x]~ 3 For the time and spatial derivatives to have the same units in £ we take 

<9o = - |r, (68) 
c at 

and so <9 M ~ This gives [/i] ~ l/y/[x] [t]. We require [m fc c 2 } ~ ~ 1 / [t] , so that 

we have [rrih] ~ [t]/[a;] 2 , and choose the electroweak symmetry breaking expectation value 
to have the same dimensions as the field h, [v] ~ l/y/[x] [t]. The La grangian density with 
these unit conventions is given by 



2 



tth tth 2 

-i . eff IL . eff IL 
l + C{ J - + cJ J — 

V v z 



1 v;A e// A e// / v 2 



Now consider the non-relativistic limit of this theory. The interaction terms will be 
determined below by matching. Consider first the theory of a free real scalar field of mass 
rrih given by 

C= l -d^d^- l -m 2 h c 2 ^ 2 . (70) 

The field tp must also be expanded in the c — >• oo limit. We remove a large energy scale 
rrih c 2 from this field with a field redefinition 

<p( x ) = e - imhc2r - x <p+(x) + e imhc2r - x <p.(x), (71) 

where r = (1,0) and <p + (x),<p-(x) correspond to the creation and annihilation components 
of the scalar field <p(x). Expanding the Lagrangian density in terms of <p±(x) we neglect 
terms multiplied by factors of 

e BimkA , (n^0,nel). (72) 
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These are terms in the Lagrangian density where some of the fields are far off shell. Their 
effect is only to modify coefficients of local operators in the effective action. 
With this substitution we find 



C = <9 y?_ d° (p + — di ip_ &" (p + + im h c (y?_ d° (f + — cp + d° y?_) . (73) 
The first term, with two time derivatives, is suppressed by 1/c 2 and is suppressed in the 



c — > oo limit. Integrating by parts the remaining kinetic terms and rescaling h± = \ / 2rrih >~p± 
gives 

£ ™ = A -( ! I + C' I+ ' (74) 

One can extend this effective Lagrangian by adding interactions, including higher order 
terms suppressed by |u|/c and v 2 /A4 2 where u is the relative velocity of the two higgs in a 
non-relativistic bound state. Scattering is described by a contact interaction which can be 
parameterized by a coupling Cnr, 

There is no cubic interaction because this necessarily involves at least one far off shell 
particle. The effect of the cubic interaction in the HET is incorporated in the coupling Cnr, 
and we will compute this in terms of the parameters of the HET below, in Sec. IVI Al 
In this effective theory, the energy and the momenta of the system are given by 

k° = l -m h \u\\ (76) 
q = m h u, (77) 

where |u| <C c is the relative velocity of the two higgs. It is advantageous to have power 
counting rules in |u| that are as manifest as possible in the Lagrangian density as demon- 
strated in 62j . We rescale so that the natural sizes of the coordinates are given by the above 
energy and momentum and define a new field H±(x) and new, dimensionless coordinates X 
and T by 

x = A a .X, t = X t T, h ± {x) = \ h H±{x). (78) 



23 



To ensure the rescaled energy and momenta are of order unity we have At = rrih\ x and 

(79) 

(80) 

(81) 
(82) 

(83) 



\ — 


1 


rrih u ' 


M — 


1 


m h u| 2 ' 


K = 


mf |u| 3 / 2 


K° = 


k° 


1 19' 

rrih u \ z 


K = 


q 


m h u 



The form of the Lagrangian density when we implement these re-scalings and introduce an 
appropriately rescaled contact coupling, Cnr = 4 ml Cnr is given by 

Cnbh = H_ (id* + ?P) H + + ^ M (i/_) 2 (i/ + ) 2 . (84) 

This form of the Lagrangian makes power counting explicit in the small parameter u/c. 
Physical quantities, such as the energy of bound states, can be equally calculated from the 
theories in Eqs. ( 1751) or ( 1841) . Which is used is a matter of convenience: the former has 
familiar dimensions while the latter has explicit power counting. 5 

A. Matching onto NRHET 

To determine the matching coefficient Cnr we take the non relativistic limit of the hh — > 
hh scattering determined in HET. We neglect the running from ml down to our matching 
scale /i 2 = m? h in this initial study, and perform the matching at tree level only. 

1. Linear Realization 



The HET contact interaction is given by 



^ = -3A 1 + 20A lC f-^A 2 ^ + i^M). (85) 



5 The c — > oo limit of NR effective field theories was studied in [63|. The reader interested in bound 
states at threshold in NRHET would also profit from an examination of the treatment of bound states at 



threshold in NN effective field theory, reviewed in 64 1 
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The Yukawa exchange feynman diagrams, shown in Fig. [5l give the amplitude 

i A^s, t, u) = (t) + A[(u) + (86) 
where s,t,u are the usual Mandelstam variables, and 

Afo) = ~K X 1- f 3 Ax + 5 A 2 £ - 4 (* + 2 ml) % - ^ f ^ V 18 A i ( 87 ) 

\ / v ^ / "-. \ ' 



+ ix + T x 

FIG. 5: Tree level hh — > hh scattering in the extended higgs theory. Time flows left to right. 
The total amplitude for hh — > hh scattering is given by 

A^ hh (s,t,u)=A^ + Ai(s,t,u). (88) 

To perform the matching we take the momenta of the higgs particles to be off-shell by a 
small residual momenta p with energy and momenta that scale as po ~ rrihU 2 and p ~ m^u. 
The momenta of the higgs are decomposed as (recall r = (1,0)) 

p = rrihr + p, k = rrihr + k, 

p' = m h r + p', k' = m h r + k'. (89) 
This gives, in the center of mass frame 

s = 4m 2 h + 4 |q| 2 , 

t = -|q| 2 (l-cos(#)), (90) 
u — — |q| 2 (1 + cos(0)) , 

with q ~ nihil. In the non relativistic limit we retain the lowest order in |u| and we have 



A L NR = A^ + Af(Aml)+2A^0) 



12A 1 + 10A 2 ^-64A 1 Cf-^ m. (91) 
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To determine the coupling C% R in the NRHET Lagrangian, Eq. (175]) . we compute the four 
point amplitude and insist that it equals Anr- Inserting four factors of y/2rrih to account 
for relativistic normalization of states, we finally arrive at 

{2m h fC L NR = C L NR = 12 A x + 10 A 2 ^ - 64 A, Cf - ^ ^ . (92) 
2. Non-Linear Realization 

For a non-linear realization we find the following for the HET contact interaction 



2 



Aq L = -X e / f + 4?±d?'. (93) 
The yukawa exchange diagrams give 



A? L {x) = 1— 

x — vn h + i e 



c/ f ( 2m\ + x 



1 2 



(94) 



2 \ v 2 

The matching is performed as in a linear realization and we find 
A N N L R = A» L + A? L (±ml)+2A? L (0) 

This gives the effective HET coupling in the non-linear realization 

{2m h fC N A = 6™ = \- 2 {\r? - K" ~ + i^f - (4"r) 4- 06) 

6 mi V / v z 



B. NRHET Bound State Energy 

To find the approximate bound state energy of the higgs, we calculate the bubble sum 
in our NRHET theory and interpret the pole in the re-summed bubble chain as the bound 
state energy of Higgsium. Note that this calculation is formally justified in the large N limit 



651 ] where the higgs sector is equivalent to and 0(4) theory [66|. The Feynman rules for the 
NRHET Lagrangian in (!75|) are shown in Fig. [61 

The bubble sum is straightforward to calculate in NRHET. The leading order term is 
directly obtained from the Feynman rules, we use the Lagrangian given by Eqs. (1751) in the 
following. The leading bubble graph is given by 

n fdk d d k 



iAl-l 0ov -{iC NR ) I {27r)d {E + k0) _ k 2 /2mh + , e - _ k 0_ k 2 /2mh + te ( 97 ) 
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FIG. 6: Feynman rules for NRHET. 

We have chosen to work in the center of mass frame, and E = + P® stands for the center 
of mass energy. Performing the first integral by residues and the remaining integrations with 
dimensional regularization, we find 



iAi- 



loop 



m h (C_ 



NR) 



4lX 



-m h E) 



1/2 



(98) 



+ 



+ 



FIG. 7: The bubble sum of graphs leading to the bound state pole in NRHET. 

The terms in the bubble sum of diagrams shown in Fig. [7] are given by the geometric 
series 

2 



iC 



NR 



1 - 



m h C NR , ^M/ 2 , ( rn h C NR mhE y/2 



An 



\ 4:71 



+ 



iC 



NR 



1 I m hCNR | 

1 " r 4tt ' 



-m h E)V* 



This result agrees with 62J, |67J and indicates a bound state with a bound state for Cnr > 
with binding energy 

1 / Ait \ 2 f 1fW \ 2 

(99) 



1 ( 4?r 



Vm h C, 



C 



NR 



There is an implicit renormalization condition introduced by dimensional regularization. 
The integral has no pole as d — > 3, so it is interesting to ask what subtraction has been 
made. This is easily understood by performing the d = 3 integration with a momentum 
cut-off |k| < A in terms of the bare coupling C% R : 



iAt loop = im h (C? 



'NR) 



A 



2tt 2 An 



(100) 
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The renormalized coupling CVr(/i) can be denned as the amplitude at a fixed energy E = —jj, 



681 ] . Then the combination 



1 1 m h A 1 (rrihu) 1 / 2 

- m h K h -f> (101) 



Cnr C° nr 2tt 2 C NR {fi) 4vr 

is renormalization group invariant. This is precisely the coupling that appears in (l98l . 

It would appear that for any positive value of Cnr we have bound states. However for 
our NR description to be self consistent we require that the binding energy of the bound 
state satisfy E b < m h , that is, 

C NR > 16?r. (102) 

1. Linear Realization 

In the case of a linear realization a heavy higgs seems necessary for the bound state 
to form, but the new physics effects may allow significantly smaller masses for the bound 
state. If we neglect the effects of new physics (A2 and Off) the bound (11021) translates into 
nth > 2.0 v. Retaining the effectsof A2 and C^ one determines a condition for the NRHET 
calculation of the bound state energy to be self consistent 



v 



16ir-4A 2 & + %fe(£ 
\J 12-40C* 



(103) 



Alternatively, for a given value of the higgs mass, say ra^ = £i>, the self-consistency 
condition implies a constraint on the coefficients of the higher dimension operators: 

12 « 2 + 4 ^ < A . " 10 « 2 < C J + \0) > + 7# (^) HO 4 ) 
Using M. = 1 TeV and the PDG value for the top quark mass, this condition simplifies to 

1.2 + 0.024A! - 0.24 £ 2 (Cj + ^Cj) > 5.1 (105) 

So, for example, for \C^ + |C?| = 1 or 5 the minimal higgs mass for a NR bound state is 
reduced by 6% or 28%, respectively. Near the limit of validity of our calculation ~ 2 rat, 
for negative values of C\ + ^C? we find that a bound state is possible for Oil) wilson 
coefficients as we illustrate in Fig 7. 
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FIG. 8: In the linear realization the allowed parameter space for NR bound state formation is 
above the line. 

2. Non-Linear Realization 

In the nonlinear realization this condition is easily satisfied even for a light higgs, mh < v. 
Recall that X e / f and c\ are both enhanced by powers of Ai/v. Taking, for example, 
m h = 120 GeV and Ai — 1 TeV, neglecting the contribution of , the NR bound state 
condition is 

5 M 2 ~ eff 2 ~e//x~ e // m h (~ eff\2 . n R „ , x / 1nfi N 

- — T (A 3 ) -2ci-'''A3 — 779 ( Cl ) >lo7r + A 4 , (106) 

where 

a;« = (^) a* , 

<«" = (jz) (107) 

Note that as M. grows larger the region that satisfies the NR bound state condition grows. 
This is due to the fact that the attractive interaction given by \ e J f rsj M/v is a relevant 
operator. We find that as grows and as Ai is larger the allowed parameter space of the 
NR bound state condition is significant, demonstrating that a bound state is likely to form 
in the non linear realization. 
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m h = 12 GeV 
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FIG. 9: In the Non-linear realization, holding fixed and set A4 = as it is 0(1) and suppressed 
by 16 7r. We vary M for the values A4 = ITeV (dotted line), A4 = 3TeV (dashed line) and 
M. = 10TeV (solid line). The region above (the upper) and below (the lower) hyperbolic curves 
satisfy NR the bound state condition. 
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FIG. 10: In the Non-linear realization, holding M fixed and set A4 = as it is 0(1) and suppressed 
by 16 7T. We vary m/j = 300 GeV (solid line), = 200 GeV (dashed line) and = 100 GeV 
(dotted line). The region above (the upper) and below (the lower) hyperbolic curves satisfy the 
NR bound state condition. 

VII. SUMMARY AND CONCLUSIONS 

If a new strong interaction is responsible for electroweak symmetry breaking but a higgs 
particle, the pseudo-goldstone boson of broken scale invariance, remains unnaturally light, 
the self-interactions of this higgs particles could be quite strong. If strong enough these 
self-interactions could bind two higgs particles. 




30 



To study these questions we formulated two different effective theories of the light, self- 
interacting higgs below the scale A4 of the new physics. In the first, the symmetry is realized 
linearly and the higgs field is described as one component of an SU(2)^ doublet, just as in 
the standard model of electroweak interactions. In the second approach the symmetry is 
realized non-linearly: the triplet of would-be goldstone bosons and the higgs field are not in 
a common multiplet. We note that operators of dimension 3 in the effective Lagrangian in 
the non-linear realization are naturally expected to be enhanced by a power of A4/v relative 
to their linear realization counterparts. 

In order to study how large these couplings need be, we have studied the case of non- 
relativistic bound states. To this end we constructed a non-relativistic higgs effective theory 
(NRHET) describing self-interacting higgs particles in the rest frame of the bound state, in 
the non-relativistic limit. 

The effects of the top quark are small but non- negligible. We estimated them by including 
the virtual top quark effects as a modification to the couplings in the NRHET. 

Our results show, perhaps not surprisingly, that in the non-linear realization it is quite 
easy to form light Higgsium, as we call the higgs-higgs bound state. For natural couplings 
in the linear realization a bound state is only likely to form for rrih ~ v. Relativistic bound 
states are possible in both the linear and nonlinear realizations. 

There are many questions that we have not addressed. The most immediate one is how to 
search for Higgsium. Assuming a light higgs is found, one could imagine strategies involving 
invariant mass distributions of higgs-pair production. A dedicated study is required to 
determine if this or other strategies are viable. Another, related question is whether the 
effects of a short lived bound state could be seen indirectly, much like would-be toponium 
affecting the line shape in top quark pair production near threshold in e + e~ collisions. 
It would also be interesting to solve the bound state equation in the more general, fully 
relativistic case. We hope to return to these problems in the future. 
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APPENDIX A: CUSTODIAL SYMMETRY AND THE S PARAMETER 



There is some confusion in the literature regarding custodial symmetry and the operator 
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iWfiB^Wipv (Al) 



M 2 

which corresponds to the S parameter. Consider the matrix representation of this operator 



691 ] where the Higgs doublet field is given by 



(A2) 



Then e<p* is also an SUl(2) doublet with components 



(A3) 



where <p = <ft + *. The Higgs bi-doublet field is given by 



1 



V2 y-<j,- 

The SUl(2) x Uy(1) gauge symmetry acts on the Higgs bi-doublet as 



(A4) 



SU L {2) : $^L0 (A5) 
U Y (1) : $^$e^ a39/2 . (A6) 

In the limit that hyper charge vanishes the Lagrangian also has the following global symmetry 

SU R {2) : $-^i*t. (A7) 

When the Higgs acquires a vacuum expectation value, both SUl(2) and SUr(2) are broken, 
however the subgroup SUl=r(2) is unbroken, ie 

L($)L t = ($). (A8) 

This is explicitly the custodial symmetry, and the corresponding transformation of the Higgs 
bi-doublet under this symmetry It is easy to see that 
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M 2 
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Tr {$>WWip V &) B» u (A9) 



is invariant under this symmetry. The Higgs bi-doublet transforms as above and the field 
strength a I Wi fJiU transforms as 

o l W Iilv ^Lc l W Iilv l). (A10) 

However, it is also easy to see that this representation of the operator vanishes by explicitly 
performing the trace; one finds 

Tr(&a I W IliV $) = 0. (All) 

The non trivial representation of the operator in terms of the bi-doublet is given by 

- Tr ($t a 1 $ a 3 ) . (A12) 

With this factor of cr 3 , required for a non trivial representation in terms of the Higgs bi- 
doublet, one finds that this operator violates custodial symmetry 



APPENDIX B: TOP QUARK OPE 

As an example of the effect of the neglected terms in in the top quark OPE, consider the 
OPE corrections to the four point function of the higgs. The amplitude is given by 



iA 4 (s,t,u) = -6N C 



( m tY 



d d k 



Tr 



+ m t ) {k' + <fi + m t ) (tf+ft + mt) (ff + j + m t) 

mi 



\ v J J (2n) d k 2 — m\ (k + a) 2 — m 2 (k + b) 2 — m 2 (k + c)'" 



We find the leading order in p 2 /m 2 — > the amplitude is given by 



iA° 4 ( Sj t,u) 



-24 N, 



c 



mtY 



d d k (mf + 6k 2 m 2 + k 4 ) 



J J (2n) d {k 2 -m 2 f 



iN c /mty /24 _ 
16 vr 2 V v J \ e 



64 + 24 log 



V 1 

2 

mi 



(Bl) 



The leading order matching gives a factor of —4Ncmf/v 4 



Consider performing the top quark OPE to higher orders. We find for the next order in 
p 2 /m 2 

iN c (m t \ 4 ( 1 



iA\{s,t,u) 



16 7T 2 



30 m 2 t 



(a 2 + b 2 + c 2 + a ■ b + 6 a ■ c + b ■ c) . (B2) 



The invariants of the external momenta a, b, c averaged over the sum of all A4 diagrams can 
be expressed in the Mandelstam variables. We find that our momenta expressed in terms of 
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these variables are 



(a 2 ) = 41ml 

(b 2 ) = 8(s + t + u), 

(c 2 ) = 4! [(s + t + u) -3m 2 ] , 

(a ■ b) = A (s + t + u), 

(a ■ c) = 8 (s + t + u — 3 m 2 ), 

(b-c) = 16 (s + t + u-3m 2 h ) . (B3) 

With these substitutions, the next order in the expansion gives 

i N c (m t \ A ( rnW (s + t + u 3 



'^•J - (=9* (3) ( : 

i N c (m t \ A ( m 2 h \ 2 



Am\ 5 



16 ir 2 V v ) \ m 2 



(B4) 



Where in the last expression we simplified with s + 1 + u = 4m 2 l . This term matches onto 
the operator 

2 h ° = ^ 2 d»hd,h, (B5) 

with a Wilson coefficient that contains contributions from the integrating out TeV scale new 
physics and the top quark. At the scale fi 2 = m\ the Wilson coefficient is 



CfV?) = ^(4C,fK) + ^^|). (B6) 



The later term in the Wilson Coefficient is an example of a term that is neglected in our 
calculation. Corrections of this form can be systematically included by taking the top quark 
OPE to next order in p 2 /m 2 . 
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